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Abstract 

While  conventional  tomography  is  associated  to  the  Radon  transform 
in  Euclidean  spaces,  electrical  impedance  tomography  or  EIT  is  associ¬ 
ated  to  the  Radon  transform  in  the  hyperbolic  plane.  We  discuss  some 
recent  work  on  network  tomography  that  can  be  associated  to  a  problem 
similar  to  EIT  on  graphs  and  indicate  how  in  some  sense  it  may  be  also 
associated  to  the  Radon  transform  on  trees. 

1  Introduction 

As  communication  networks  have  become  an  essential  part  of  everyday  life,  disruptions 
may  have  very  serious  consequences.  Thus,  the  need  to  prevent  or,  at  least,  detect  them 
early  on,  has  become  very  important.  In  order  to  do  that  we  discuss  two  models  of 
the  problem,  one  based  on  weighted  graphs  and  the  second  based  on  trees.  The  first 
one  is  the  discrete  equivalent  of  the  inverse  conductivity  problem,  that  is,  of  Electrical 
Impedance  Tomography.  The  second  model  was  mentioned  recently  by  E.  Jonckheere 
and  his  collaborators  [29], 

In  both  cases,  the  data  we  collect  are  obtained  by  monitoring  traffic  only  at  distin¬ 
guished  subsets  of  the  network.  We  think  about  this  subset  as  being  the  periphery  of  the 
network. 

2  The  weighted  graph  model 

In  this  case  we  model  our  network  in  the  following  way.  We  have  a  collection  of  nodes 
and  edges  between  the  nodes  in  a  finite  planar  connected  graph  G.  We  denote  by  V  the 
set  of  nodes  of  G  and  by  E  the  set  of  edges  of  G.  Usually,  the  graph  G  is  denoted  by 
G(E,  V ).  A  particular  subset  of  this  graph  G  is  denoted  by  dG  and  called  the  boundary 

1  The  present  is  an  extended  version  of  the  lecture  given  by  Prof.  Berenstein  at  the  special  session  on 
tomography,  2004  AMS  meeting.  Rider  University.  The  authors  acknowledge  partial 
support  on  this  research  from  grants  ARO-DAAD-1901 10494  and  NSF-DMS-0400698. 

3  Department  of  Mathematics,  University  of  Maryland,  College  Park,  Maryland,  20742,  USA 
3  Institute  for  Systems  Research,  University  orMaryland,  College  Park,  Maryland,  20742,  USA. 


of  G.  In  our  context  these  are  the  nodes  accessible  to  whoever  is  trying  to  monitor 
the  traffic  in  G.  The  boundary  edges  are  those  links  whose  two  endpoints  are  in  dG. 
We  assume  that  G  remains  connected  even  if  we  remove  the  boundary  edges.  For  our 
present  purposes,  the  boundary  edges  play  no  role,  thus  we  may  as  well  assume  that 
there  are  none.  We  also  assume  that  dG  is  not  empty. 

Furthermore,  we  assume  that  to  every  edge  in  E  we  have  an  associated  non-negative 
number  oj(x,  y)  which  corresponds  to  the  traffic  between  the  endpoints  x  and  y  of  the 
edge.  Note  that  this  is  a  static  model  and  we  are  really  thinking  that  the  graph  is  a  planar 
graph,  although  this  is  not  used  anywhere  in  the  reasoning.  We  define  the  degree  dux 
of  a  node  x  in  the  weighted  graph  G  with  weight  u>  by 

dux  =  ^2  u(x,y) 
yev 

the  Laplacian  operator  corresponding  to  this  weight  10  is  defined  by 

AJW  =  t/(y)  “  X  G  V 

itv  d“x 

A  graph  S  =  S( V' ,  E')  is  said  to  be  a  subgraph  of  G(E,  V)  if  Vt  CV  and  Ef  C  E. 
In  this  case,  we  call  G  a  host  graph  of  S.  The  integration  of  a  function  /  :  G  — >  R  on  a 
graph  G  =  G(V,  E)  is  defined  by 

/  /  =  V]  f{x)dulx  or  simply  /  fdu 
JG  xev  Jg 

For  a  subgraph  S'  of  a  graph  G  =  G(  V.  E)  the  (node)  boundary  OS  of  S  is  defined  to 
be  the  set  of  all  nodes  z  €  V  not  in  S  but  adjacent  to  some  node  in  S,  i.e., 

dS  =  {z  €  V  |  z  ~  y  for  some  y  G  S} 

O 

and  the  inner  boundary  dS  by 

dS  =  {z  e  S  |  y  ~  z  for  some  y  G  dS} 

where  z  ~  y  means  that  the  two  nodes  2  and  y  are  connected  by  an  edge  in  E.  Also, 
by  S  we  denote  a  graph  whose  nodes  and  edges  are  in  S  U  OS.  The  (outward)  normal 
derivative  -§£~{Z)  at  z  G  dS  is  defined  to  be 

° yes  “ 
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where  d'^z  =  £  w(z,y) 
yes 

In  this  model,  there  are  two  kinds  of  disruptions  of  traffic  data  that  could  arise. 
In  one  of  them,  disruptions  occurs  when  an  edge  “ceases”  to  exist,  in  this  case  the 
“topology”  of  the  graph  has  changed,  and  we  refer  to  the  important  work  of  Fan  Chung 
and  her  collaborators  which  offers  crucial  insights  into  this  question.  (See,  for  instance 
[16],  [17]  and  [18].).  In  the  other,  the  weights  change  because  of  “increase”  of  traffic, 
that  is,  the  network  configuration  remains  the  same  but  the  weights  have  either  increased 
or  remained  the  same.  In  this  second  situation,  we  can  appeal  to  the  following  theorem 

Theorem  1  [11]  Let  u>i  and  w2  be  weights  with  ui\  <  u>2  on  S  x  S,  G  a  graph  and 

fi,  f 2  '■  S  — >  R  be  functions  satisfying  that  for  j  =  1,2, 

f  ^ujj{x)  =0,  x  G  S 
)  Jg_(*)  =  $(*),  zeds 

1  Is  fjdwj  =  K 

for  any  given  function  $  :  dS  — >  R  with  f9S  $  =  0,  and  for  a  suitably  chosen  number 
K  >  0.  If  we  assume  that 


(i)wi(z,y)  =  0J2 {z,  y)  on  dS  x  dS 
(O)  Alas  =  /2|as> 

then  we  have 

fi  =  h  on  S 

and 


u>  1  =  u>2  on  S  x  S 

whenever  f^x)  ±  fx{y)  and  f2(x)  ±  f2(y)  ■ 


We  conclude  that  the  data  distinguishes  the  two  cases.  That  is,  we  can  decide 
whether  there  is  an  increase  of  traffic  somewhere  in  the  network  or  not.  While  this  is 
only  a  uniqueness  theorem,  nevertheless,  we  can  effectively  compute  the  actual  weights 
from  the  knowledge  of  the  Dirichlet  data  for  convenient  choices  of  the  input  Neumann 
data  in  a  way  similar  to  that  done  in  [21]  and  [23]  for  lattices.  Similarly,  the  Green 
function  of  this  Neumann  boundary  value  problem  can  be  represented  by  an  explicit 
matrix. 

What  we  want  to  discuss  now  is  the  relationship  between  the  above  results  to  the 
problem  of  understanding  a  large  network  like  the  internet.  One  way  to  make  more 
concrete  this  problem  was  discussed  by  T.  Munzner  in  [32]  and  [33]  on  visualizing  the 
internet.  It  implies  that  the  natural  domain  might  be  a  hyperbolic  space  of  dimension 
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higher  than  2.  One  can  see  that  Munzner’s  suggestion  leads  to  a  question  closely  re¬ 
sembling  EIT,  and  it  is  natural  to  consider  it  a  problem  in  hyperbolic  tomography  [7], 
[8].  On  the  other  hand,  we  have  just  obtained  a  significant  result  on  the  inversion  of  the 
Neumann-Dirichlet  problem  by  studying  it  directly  on  “weighted”  graphs  [11],  Sim¬ 
ilarly,  the  Radon  transform  in  the  hyperbolic  plane  has  been  studied  in  [7],  [8],  and 
[27], 

In  addition,  in  a  recent  lecture,  E.  Jonckheere  [29]  indicated  that  at  least  locally  in¬ 
ternet  traffic  could  see  modelled  as  being  part  of  a  tree  and  therefore  it  can  be  visualized 
using  2-dimensional  hyperbolic  geometry.  As  a  consequence,  a  different  way  to  study 
locally  this  kind  of  networks  can  be  done  using  the  Radon  transform  on  trees.  As  it 
turns  out,  inversion  formula  for  the  Radon  transform  on  trees  is  already  known  and  it 
can  be  found  in  [9], 

For  the  sake  of  completeness,  we  will  describe  here  a  simplified  version  of  the 
Radon  transform  on  trees  and  its  inversion  formula.  As  explained  below,  this  seems  to 
be  enough  to  deal  with  the  network  problems  we  are  interested  in. 

3  The  Radon  transform  on  homogeneous  trees 

Let  us  now  remind  the  reader  what  do  we  mean  by  a  tree  T.  A  tree  T  is  a  finite  or 
countable  collection  Uof  vertices  {vj,  j  =  0, 1, ....}  and  a  collection  E  of  edges  ejk  = 
(■ Vj,Vk ),  in  other  words,  pairs  of  vertices.  We  orient  the  edge  ejk  by  thinking  that  Vj 
is  the  first  node  and  17  the  second  node.  We  always  include  the  edges  e\ y  in  this 
collection,  which  have  the  reverse  orientation.  Given  two  vertices  u  and  r.we  say  they 
are  neighbors  if  (it,  v)  is  an  edge  and  write  u  ^  v  in  this  case.  A  geodesic  7  from  uq  to 

ui  is  a  collection  uo  ,  '«i  , . ,  ui- 1 ,  ui  of  pairwise  distinct  vertices  such  that  uq  ^  u\  , 

u±  ^  U2  , _ ,ui—\  ^  ui  ■  It  turns  out  that  uq  ^  Ui  then  we  consider  the  closed 

geodesic  path  7  by  adding  the  edge  (iq  ,uq)  to  7.  Unless  explicitly  mentioned,  our 
geodesics  will  not  be  closed.  To  simplify  the  notation,  for  any  geodesic  7  =  u0 
Mi  ^  Mi  ^  U2  ^  ••••  ^  M;_i  ^  ui  open  or  closed,  we  denote  by  —7  the  geodesic 
with  the  opposite  orientation  ,  i.e.,  -7  =  u  1  ^  m;_i  ^  ....  ^  uq  .  The  collection  of 
all  (open)  geodesics  is  denoted  by  T.  If  T  is  infinite,  then  a  complex  valued  function 
f&L\T)  if. 

i/wi  < 00 

vev 

the  Radon  transform  li  of  a  function  f  £  l?  (1")  is  simply  the  bounded  function  Rf  on 
r  defined  by 

Rfil)  =  Y  /(v) 

Given  a  node  v  we  denote  by  v  ( u )  the  number  of  edges  that  contain  v  as  an  endpoint. 
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This  number  is  sometimes  called  the  degree  of  the  node.  We  will  assume  throughout 
that  we  always  have  v(v)  >  3  to  ensure  that  the  Radon  transform  in  injective.  (In 
our  applications  this  is  only  needed  for  nodes  v  that  lie  in  supp(/).  In  the  terminology 
of  [9]  we  are  assuming  there  are  neither  black  holes  nor  flat  points  in  T.  Under  these 
conditions,  the  Radon  transform  in  a  tree  is  invertible.  In  fact,  the  explicit  inversion 
formula  resembles  that  of  the  inversion  for  the  Radon  transform  in  the  Euclidean  plane 
[10],  [12],  [13],  and  [27].  Unfortunately,  even  in  this  case,  we  need  to  introduce  a 
significant  amount  of  auxiliary  notation.  For  the  purpose  of  illustration  we  describe  the 
inversion  formula  here  only  for  the  case  of  homogeneous  trees. 

4  Inversion  of  the  Radon  transform  in  homogeneous 
trees 

Consider  a  homogeneous  tree  T  in  which  each  vertex  touches  q  +  1  edges  with  q  >  2. 
If  n  is  a  nonnegative  integer,  let  V(n  the  number  of  vertices  of  T  at  distance  n  from  a 
fixed  vertex  of  T.  It  follows  that 


f  1  if  n  =  0 

\  (q  +  l)gn_1  ifn  >  1 

We  give  the  following  definitions.  Let  v,w  two  vertices  in  T  that  are  connected  by 
a  path  (v  =  vo,  =  w),  then  the  distance  between  v  and  w  is  the  nonnegative 

integer  |u,  w\  =  m.  Also,  for  /  £  L1(T),  let  pn  the  average  operator  defined  by 

hnf(v)  =  —  /M>  for^eT 

\vM=n 

It  can  be  seen  that  fin  is  basically  a  convolution  with  radial  kernel 


hn(v,w)  = 


v(« ) 

0  if 


|u,  w\  =  n 
|u,  w |  ^  n 


Let  P  =  q/(2(q+  1))  and  R*  be  the  dual  Radon  transform  defined  for  $  £  L°°(T) 


by 

/?,*<!> (v)  =  J  for  each  vertex  v  £  T, 

with  respect  to  a  suitable  family  {pv  :  v  £  T}  of  measures  on  I’  where  I is  the  set  of 
all  of  the  geodesics  containing  the  vertex  v. 

In  order  to  obtain  the  inversion  of  R  we  observe  that  R*R  acts  as  a  convolution 

OO 

operator  given  by  the  radial  kernel  h  =  (3 ho  +  ^  2 (3hv. 

n— 1 
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Proposition  2  The  identity 


R*R  =  /3 fio  +  onL\T), 

n—1 

holds  in  L 1(T),  where  the  series  is  absolutely  convergent  in  the  convolution  operator 
norm  on  L2(T),  thus  providing  a  bounded  extension  of  R* R  to  L2(T). 

Theorem  3  The  unique  bounded  extension  to  L2(T )  of  the  operator  R*R  is  invertible 
on  L2(T),  and  its  inverse  is  the  operator 

OO 

To  +  y~!(-l)w2  pn 

n= 1 

which  acts  as  the  convolution  with  the  radial  kernel  [ho  +  ^2  (— l)n2 hn\.  As 

n=  1 

before,  this  series  converges  absolutely  in  the  convolution  operator  norm  on  L2(T );  in 
particular,  E  is  bounded. 

Corollary  4  The  Radon  transform  R  :  Ll(T )  — >  Loc(Tj  is  inverted  by 

ER*Rf  =  f. 


2(g  +  l)3 

q{q~  l)2 
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